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Resume 

La technique transaurale est une methode qui permet de restituer un 
environnement sonore spatialise avec des enceintes en utilisant un signal binaural. 
Elle met en jeu differents filtrages destines a annuler les chemins croises tout en 
preservant le timbre de I'enregistrement. Ce rapport de stage presente mon travail 
d'optimisation du procede et sa realisation pratique au sein de la societe GENESIS. 


Transaural acoustic is a method to reproduce a 3 dimensional sound with 
several speakers using a binaural recording. It involves successive filtering in order 
to cancel crosstalk without disturbing the recording's timber. This report deals with 
my work on the optimization and the realization of transaural during my internship 
in GENESIS. 


Mots cles : 


Binaural, Basses frequences, Conditionnement de matrice, Egalisation, HRTF, Inversion, 
Localisation, Matrice de transfert electro-acoustique, Optimal Source Distribution, 
Regularisation, Stereo-dipole, Timbre, Transaural. 
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Glossaire 

Conditionnement de matrice : Mesure de la difficulty du calcul numerique du probleme. 

Enregistrement Binaural : Methode cherchant a reproduire la perception sonore naturelle 
humaine par restitution au casque. 

HRTF (Head related transfer function) : Donnees caracterisant les transformations 
apportees a I'onde sonore par le corps d'un auditeur. 

Inversion de Kirkeby : Methode d'inversion de filtre d'egalisation utilisant la regularisation. 

Matrice de transfert electro-acoustique : Matrice qui caracterise tout la chaine de 
propagation sonore des enceintes aux oreilles de I'auditeur. 

Regularisation de Tychonoff : Methode utilisee dans la resolution de problemes qui ne sont 
pas bien poses. 

Technique transaural : Technique qui permet de recreer un environnement sonore 
tridimensionnel avec des haut-parleurs a partir d'un signal binaural. 
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1 Introduction 


Depuis quelques annees un engouement pour les technologies de restitutions 
tridimensionnelles est ne dans de multiples domaines : 

- Les simulateurs d'environnement sonore pour I'automobile et 
I'aeronautique. 

- Les jeux video. 

- Le cinema. 

- La musique. 

A ce jour plusieurs pistes de recherches sont en cours et on a vu naitre 
diverses techniques telles que I'holophonie, I'ambisonic, la WFS, les restitutions 
binaurales et transaurales. L'objet du stage concerne cette derniere qui consiste a 
reproduire un champ de pression determine au niveau des deux conduits auditifs de 
I'auditeur avec des haut-parleurs. 


La technique est eprouvee et utilisee par GENESIS depuis quelques annees 
mais certains points restent optimisables. GENESIS utilise un systeme de 4 haut- 
parleurs + 1 caisson de basse, mais les basses frequences etant peu localisables et 
le procede de transauralisation delicat dans cette gamme de frequence, la technique 
actuelle n'effectue pas de transauralisation en dessous de 70 Hz. Lors de 
Legalisation finale en timbre des systemes en basses frequences, ceci peut poser 
des problemes pour la maitrise des deux champs acoustiques distincts au niveau 
des oreilles. 


L'objectif du stage est ainsi d'optimiser les systemes transauraux et de mieux 
maitriser les basses frequences. 


Apres une breve presentation des enjeux industriels du syteme transaural, 
nous nous attacherons dans un premier temps a la presentation de la technique en 
mettant en exergue les points defaillants. Dans un deuxieme temps nous 
presenterons les solutions apportees aux problemes precedents. 
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2 Presentation de GENESIS et des enjeux du transaural 

2.1 La societe GENESIS 

GENESIS est une societe anonyme fondee en 1999 au capital de 123 000 € 
specialise dans I'innovation en acoustique. La societe reunit 10 personnes sur le 
Technopole Europeen de I'Arbois base a Aix-en-Provence. 


GENESIS est un acteur implique dans la recherche puisqu'elle investit 30% de 
son chiffre d'affaire dans la recherche et le developpement. De plus elle possede un 
partenariat avec le Laboratoire de mecanique et d'acoustique de Marseille. 


GENESIS travaille avec de nombreux industriels dans les secteurs automobiles 
et aeronautiques comme Renault et Airbus. La societe propose des expertises et des 
produits logiciels pour contribuer a ^amelioration de la qualite sonore dans les 
moyens de transports, les environnements quotidiens (station de metro, nuisance 
des avions, bruist des chantiers...), des simulateurs interactifs et 3D (simulateurs de 
vol, de conduite, de sous-marin, des systemes audio de realite virtuelle) ... 


2.2 Le transaural 

Les industriels de I'automobile et de I'aeronautique s'equipe de salle de reecoute afin 
de faire des etudes comparatives d'environnement sonore facilement. Cette technique est 
preferee a I'ecoute binaurale car elle permet une meilleure externalisation des sources 
sonores sans pour autant avoir a par exemple conduire successivement 2 voitures pour 
comparer le son produit dans I'habitacle. Elle permet done a la fois d'etre un outil de 
developpement pour les ingenieurs qui ont acces a une comparaison facile mais aussi d'etre 
un argument convaincant les decideurs par ecoute comparee. 


GENESIS ne connaTt a ce jour aucune concurrence sur le territoire frangais, mais les 
exigences en termes de spatialisation des industrielles poussent la societe a optimiser le 
procede pour gagner en efficacite. 


Afin de satisfaire les objectifs de la mission de mon stage j'ai etudie et realise un 
transaural par la technique OSD (Optimal source distribution) qui permet d'avoir une 
meilleur spatialisation. J'ai egalement mis en place une technique d'inversion par la 
methode de Kirkeby qui permet d'avoir un resultat en timbre superieur a ce qui etait utilise 
auparavant. 
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3 Theorie des systemes transauraux 

3.1 La perception spatiale d'une scene sonore et la technique 
binaurale 

3.1.1 La perception spatiale d'une scene sonore 

La perception naturelle de I'espace sonore passe par le positionnement des 
sources qui nous entourent. Localiser une source revient a determiner sa direction, 
done son azimut et son elevation, ainsi que la distance a laquelle elle se trouve (ou 
profondeur). Ces elements sont definit sur la figure suivante : 



La localisation naturelle 


Figure 1 : La perception spatiale d’une scene sonore 


Lorsqu'une onde acoustique issue d'une source S parvient aux oreilles de 
I'auditeur et que cette source est situee par exemple sur la droite, on peut definir 
par Sd I'onde arrivant directement a I'oreille droite, et par Sg I'onde arrivant a 
I'oreille gauche de maniere indirecte, appelee alors " onde de contournement 
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l a localisation nature-lie laterale 

Figure 2 : L'onde de contournement 


La bonne localisation des sources acoustiques dans un environnement sonore 
passe par ^interpretation des differences entre ces 2 signaux Sd et Sg (1). Les 
differents indices analyses par le cerveau pour la localisation des sources dans 
I'espace sont presentes ci-apres : 


• La difference interaurale de temps (ITD) : 

La presence d'une difference de marche entre les trajets acoustiques de la 
source a chacune des deux oreilles implique une difference de temps d'arrivee pour 
le son a chacune des oreilles. On peut assimiler la tete a une sphere de rayon r et 
ainsi definir la difference de temps d'arrivee des signaux provenant d'une source 
sonore suffisamment eloignee en fonction de son angle d'incidence : 



Figure 3 : Difference de marche de l'onde contournement par rapport a l'onde directe 
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• La difference interaurale de phase : 

Elle est surtout significative en basses frequences (en hautes frequences les 
periodes sont trop petites pour que les dephasages soient interpretables). La figure 
4 montre le lien entre difference de temps d'arrivee et dephasage : 



Figure 4 : Illustration de la difference de temps d'arrivee (ITD) et de la difference de phase entre les 2 

oreilles 


• La difference interaurale de niveau (ILD) : 

Elle caracterise la difference de niveau sonore entre les deux oreilles. De par 
sa morphologie, la tete agit en effet comme un obstacle masquant et diffractant. 
Les differences interaurales de niveau dependent de la frequence, comme le montre 
la figure 5: 
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DifTertiKC IdB) 



Figure 5 : Evolution de la difference interaurale d'intensite en fonction de la position de la source sonore 
pour differentes frequences 


• Les HRTF (Head Related Transfer Functions) : 

La connaissance des 3 indices mentionnes ci-dessus n'est pas suffisante pour 
permettre une localisation parfaite des sources. En effet il demeure un cone de 
confusion puisqu'une source situee a I'avant et une source situee a I'arriere avec le 
meme angle par rapport a I'axe interaural ne peuvent etre discriminees : 



Figure 6 : Illustration du cone de confusion 


Yannick Busawon 
Avril 2013 
Page 12 / 45 


Iqenesis 

M “all dimensions of Sound 


GENESIS 

Restitution des Basses Frequences et optimisation 
dans les systemes Transauraux 3D 


Ceci conduit a ('introduction des HRTFs (Head Related Transfer Functions). En 
effet compte tenu de la geometrie des pavilions auditifs et de notre tete, le spectre 
du signal regu au tympan varie suivant la direction d'incidence. L'ensemble des 
filtres crees par la tete en fonction des differentes incidences constitue un indice 
supplementaire analyse par le cerveau humain pour localiser une source dans 
I'espace et permet notamment la discrimination avant arriere (2) (3). II faut 
cependant noter que les HRTFs sont des fonctions personnels puisqu'elles 
dependent de la morphologie de chacun. 

Pour illustration, la figure 7 montre un exemple de base de donnees de HRTF 
mesurees au M.I.T sur une tete artificielle: 


HRTFs 

(oreille gauche/plan horizontal) 



Figure 7 : HRTF mesurees au MIT par pas de 5° dans le plan horizontal 


Pour realiser cette figure 7 on mesure I'intensite acoustique regue par la tete 
artificielle sur toute la gamme de frequence en faisant tourner la source acoustique 
par pas de 5° dans le plan horizontal. Les couleurs represented ici I'intensite recue 
pour une frequence et un azirmut de la source donne. 


Les indices de perception dans le plan vertical sont moins faciles a mettre en 
evidence. On peut cependant noter que les experiences de Blauert donnent le jour a 
la notion de « bandes directives ». II a pu mettre en evidence le lien entre le 
contenu spectral d'un signal et la direction de I'espace sonore a laquelle il est 
associe. Par exemple, des sons possedant beaucoup d'energie autour de 8000 Hz 
donnent la sensation de provenir du haut. Ainsi Blauert suggere de placer les indices 
spectraux au coeur de la comprehension des mecanismes de perception de I'espace 
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sonore. Cette remarque conduit une nouvelle fois a la description du modele de 
perception a partir des HRTFs (Head Related Transfer Functions). 


La figure 8 presente les resultats d'une experience de Blauert concernant la 
probabilite de localisation d'une source selon 3 directions (dessus, devant, derriere) 
en fonction de la frequence : 



Figure 8 : Les bandes directives de Blauert : Probabilite de localisation selon une direction donnee en 
fonction de la frequence 


On constate par exemple qu'une source avec un contenu spectral fort autour 
de 1 kHz aura de fortes probabilites d'etre localisee derriere, alors qu'une source 
avec un contenu spectral fort autour de 3 kHz sera plutot localisee devant. 


3.1.2 La technique binaurale 

La technique binaurale est la methode de restitution 3D qui s'approche le plus 
de I'ecoute reelle puisqu'elle s'appuie sur les caracteristiques psychophysiologiques 
de I'audition humaine pour localiser les sons. Le principe est de recreer le champ de 
pression acoustique d'une scene sonore reelle au niveau des tympans d'un auditeur 
afin de donner toutes les informations necessaires pour la reconstruction du 
paysage sonore. 


En pratique on realise un enregistrement en utilisant une tete artificielle (qui 
modelise un individu moyen) dans laquelle on a place deux micros dans les conduits 
auditifs (on peut aussi faire cet enregistrement directement au niveau des tympans 
d'un auditeur, mais alors le signal sera optimise pour cet auditeur). La restitution 
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binaurale consiste a reecouter directement I'enregistrement au casque, avec une 
correction eventuelle pour eliminer la coloration induite par les microphones 
d'enregistrement et le casque de restitution. 



Figure 9 : Schema de principe de la technique binaurale 


Des techniques de traitement du signal permettent egalement de recreer 
artificiellement un signal binaural a partir d'une prise de son monophonique et la 
connaissance des HRTF. II suffit en effet de filtrer le signal monophonique par les 
HRTFs gauche et droite relatives a I'azimut et I'elevation que I'on souhaite obtenir 
pour synthetiser un enregistrement binaural qui aura toutes les informations 
spatiales d'une scene sonore veritable (4). 


3.2 Principe de la transauralisation 

La technique transaurale derive de la technique binaurale. Le principe est 
toujours de restituer un champ de pression determine aux oreilles de I'auditeur, 
mais la restitution se fait cette fois sur hauts parleurs (cf figure 10). 

Le passage sur haut-parleurs n'est toutefois pas trivial. En effet dans le cas de 
I'ecoute binaurale chaque ecouteur est dedie a la synthese du champ acoustique au 
niveau d'une seule oreille. Alors que dans I'ecoute transaurale il existe des chemins 
croises entre les haut-parleurs et les oreilles qui viennent perturber la creation du 
champ acoustique. L'objet de cette technique est done d'annuler ces chemins 
croises par des precedes de filtrage, comme il sera explique plus en detail dans les 
sections qui suivent. 
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3.2.1 La configuration stereo-dipole 

Considerons dans un premier temps la restitution par 2 haut-parleurs pour un 
auditeur : la configuration stereo-dipole (5). Comme nous I'avons evoque, la 
difficulty principale de cette technique est I'annulation des chemins croises afin de 
retrouver au niveau des oreilles le signal binaural d'origine. Afin de mieux 
comprendre utilisons le schema suivant : 


U1 U2 



Figure 1 1 : Mise en evidence des chemins croises 
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Si on envoie directement les signaux U1 et U2 aux hauts parleurs VI et V 2 on 
s'apergoit que les signaux qui arrivent aux oreilles (W1 et W2) sont la somme des 
contributions de chaque haut-parleur : 


Soit : 


W ± = CM + C 2 V 2 
W 2 = CM + CM 



( 1 ) 


( 2 ) 


II est done necessaire d'annuler les chemins C2. Pour se faire, nous utilisons le 
haut-parleur de droite pour annuler le chemin croise du haut-parleur de gauche et 
inversement. On obtient alors le schema de principe suivant : 


U1 U2 



Figure 12 : Schema de principe de la technique transaurale 
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Nous introduisons des fonctions de transfert H qui bien choisies permettent 
I'annulation des chemins croises (6). Le signal envoye au haut-parleur est alors une 
combinaison lineaire des deux voies de I'enregistrement binaural. 


En introduisant les filtres H on a : 



On obtient done I'egalite matricielle suivante : 


( W A = ( c ± c 2 \(H ± h 2 \/Ua 

\w 2 ) \c 2 cJ\h 2 hJ \u 2 ) 


(3) 


(4) 


Notre but etant de recreer le signal binaural au niveau de chaque oreille, nous 
cherchons a avoir : 



Cela implique I'egalite matricielle : 

(H, H 2 \ = (C 1 CA- 1 

\h 2 hJ \c 2 cj 


(5) 


( 6 ) 


Au final I'annulation des chemins croises que nous cherchons a obtenir dans la 
methode transaurale se resume a un probleme d'inversion de la matrice de transfert 
electroacoustique entre les haut-parleurs et les oreilles de I'auditeur. Ceci pose alors 
le probleme du contenu de cette matrice, et de son conditionnement. 

3.2.2 La matrice de transfert electroacoustique 

La matrice de transfert electroacoustique contient toutes les informations 
relatives a la chaine audio (convertisseurs, amplificateurs, haut-parleurs...) et a la 
propagation des ondes acoustiques entre les haut-parleurs et les oreilles dans la 
salle consideree. On peut done etablir cette matrice a partir de la reponse 
impulsionnelle du systeme que Ton vient mesurer a I'aide d'un sinus glissant 
logarithmique ou d'un bruit rose. Cependant la mesure directe de cette matrice pose 
des problemes dans la realisation d'un bon systeme transaural. En effet les 
systemes transauraux developpes par GENESIS sont generalement installes dans 
des salles d'ecoutes qui ne sont pas anechoiques, et les reponses impulsionnelles 
mesurees presentent souvent des echos. Apres inversion nous obtenons des pre- 
echos qui viennent perturber le fonctionnement des filtres transauraux. Ce 
phenomene peut etre en partie attenue par I'introduction d'un lissage complexe 
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effectue sur les reponses impulsionnelles, mais GENESIS a choisi au final d'autres 
solutions alternatives. 


Ainsi GENESIS a pris le parti d'utiliser soit des HRTFs generiques mesurees en 
chambre anechoique, soit des HRTF issues d'un modele de tete spherique ou d'un 
modele ellipsoidal. Le choix d'une technique par rapport a I'autre depend des 
attentes formulees par le client. L'utilisation des HRTFs mesurees rend le transaural 
tres precis en localisation mais detimbre assez consequemment le signal. Le modele 
de tete spherique quant a lui detimbre moins, est plus facile a inverser et apporte 
done de la robustesse mais la localisation est moins precise. Pour finir le modele de 
tete ellipsoidale etant plus proche de la geometrie de la tete donne de meilleurs 
resultats. Cependant ce modele est analytique et la base de donnee que nous 
utilisons ne couvre le spectre que jusqu'a 5000Hz. II n'est done pas possible de 
I'utiliser sur tout le spectre. C'est pourquoi le modele utilise pendant mon stage a 
ete un modele spherico-ellipsoidal, dont I'experience a montre en interne qu'il etait 
le plus satisfaisant a I'ecoute. 


3.2.3 Le conditionnement de la matrice de transfert electroacoustique 

Le decodeur transaural etant par definition I'inverse de la matrice de transfert 
electroacoustique, le conditionnement de cette derniere est une grandeur qui par 
definition va juger de la stabilite de son inversion, et done directement de la qualite 
de notre systeme. En effet Nelson a montre que lorsque la matrice est bien 
conditionnee, cela permet a I'auditeur de faire des petits mouvements de tetes sans 
que la perception soit alteree (7). 


Definissons dans un premier temps le conditionnement de la matrice C. On se 
place dans le formalisme du stereo-dipole : 


u(k) 


^ H(k) 


vCk) 


3 C(k) 


w(k) 


Figure 13 : Representation matricielle du stereo-dipole 


Id, u represente le signal binaural en entree du systeme. H correspond au 
filtre transaural qui nous permet d'annuler les chemins croises et d'egaliser la 
reponse. v est le signal envoye aux haut-parleurs. C est la matrice de transfert 
acoustique decrite precedemment. Et w est le signal regu aux oreilles de I'auditeur. 
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Dans le cas ideal il faut : 


H = C _1 (7) 

On obtient alors v par : 

C.V = U (8) 

Cependant la matrice de transfert electro-acoustique C est issue de la mesure 
des fonctions de transfert entre les transducteurs et les deux oreilles de I'auditeur. 
La matrice C est done sujette aux bruits de mesure et cela va entrainer des 
variations sur v : 


(C + SC ). ( v + Sv) = u 


On montre alors que : 


II Sv || 

II v + Sv || 


<11 c ii. ii c - 1 


II SC || 
II C II 


(9) 

( 10 ) 


Le nombre j(C) =|| C ||. || C -1 || est le conditionnement de la matrice C, il 
se comporte comme une borne de I'erreur de la solution v qui est generee. /(C) est 
toujours superieur ou egal a 1. Si le conditionnement est egal a un le probleme est 
bien conditionne et I'inversion donne de bon resultats. Plus le conditionnement est 
important, plus I'erreur commise sur la mesure de C est susceptible de donner un 
vecteur v tres eloigne de la solution optimale. Un resultat similaire peut etre obtenu 
pour le probleme plus general a M haut-parleurs de restitution. II faut alors 
symetriser le probleme, e'est a dire trouver v qui verifie I'equation normale : 


C*. C.v = C*.u 

(* operateur de transconjugaison) 


(ID 


Dans ce cas le conditionnement peut s'exprimer en fonction des valeurs 
singulieres Oi de C : 


, ®max 

*(C) = (12) 

O min 

Nelson realise un calcul theorique en configuration stereo-dipole qui permet 
d'obtenir le conditionnement de la matrice de transfert electro-acoustique en 
fonction de I'ecart angulaire entre les haut-parleurs et la frequence (figure 14). 
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Figure 14 : Conditionnement d'un stereo-dipole en fonction de I'ecart angulaire et de la frequence de 
fonctionnement 


Les zones claires sont les zones ou la matrice est bien conditionnee. On 
s'apergoit qu'il n'est pas possible d'avoir un bon conditionnement sur tout le spectre 
avec un seul stereo-dipole quel que soit son ecartement angulaire. 


Nelson demontre theoriquement que lorsque le conditionnement est proche de 
1, le systeme possede une bonne immunite aux petites variations de position des 
points de restitution. Done si le systeme est bien conditionne, il autorise de petits 
mouvements de la tete de I'auditeur sans que la restitution trois dimensions soit 
penalisee. 

Pour realiser la figure 14 Nelson calcule analytiquement le conditionnement 
/(C) de la matrice de transfer! C pour differentes orientations de haut-parleurs. II 
s'appuie sur une representation en ondes planes du champ de pression en eliminant 
I'influence de la tete entre les deux points d'ecoute. 
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• Calcul du conditionnement pour un stereo-dipole 

Soit une source monopolaire harmonique dont les dimensions sont petites 
devant la longueur d'onde. Si q est le debit de cette source alors la pression 
acoustique en un point est donnee par la relation suivante qui est la solution de 
I'equation d'Helmholtz : 


P(M) = 


jo)p 0 q e 


-jkr 


An ' r 

Avec p 0 masse volumique, co pulsation de I'onde, k = w / c nombre d'onde, r 
distance entre la source de debit et le point M considere. 


( 13 ) 


Nelson utilise ce resultat lorsque deux sources monopolaires generent 
simultanement le champ acoustique. II cherche alors a exprimer la matrice de 
transfert electroacoustique C entre ces deux sources et deux points particuliers de 
I'espace qui seraient I'entree des conduits auditif de I'auditeur. Voici la 
representation du probleme : 



Figure 15 : Moderation du systeme a 2 monopoles et 2 points d'ecoute 

Le probleme est parfaitement symetrique. En posant g = et Al = l 2 — 

l lr nous pouvons exprimer la matrice de transfert electroacoustique C de la maniere 
suivante : 


\c n 

C 12] 

_ P 0 e Jktl 

1 ge~ ikM 

[ C 2 i 

C22I 

An' 

ge~i kM 1 


( 14 ) 
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On peut alors calculer numeriquement le conditionnement du probleme en 
fonction de la frequence et du parametre d'ecartement 0 des deux sources 
acoustiques pour obtenir la figure 14. 


• Comparaison du calcul theorique au calcul avec HRTF 

Le calcul de Nelson est theorique et ne prend pas en consideration les 
deformations du champ acoustique induites par la tete. Un calcul de 
conditionnement du probleme a ete refait en utilisant la base de donnees HRTFs du 
Massachussetts Institut of Technologies par Fabrice Gagneux a GENESIS. II obtient 
alors la representation suivante: 



Figure 16 : Conditionnement reel du probleme diversion 


Bien que Nelson ne prenne pas en consideration I'influence de la tete dans le 
calcul theorique, on peut observer qu'il obtient tout de meme une approximation 
remarquable de la realite en particulier aux basses frequences (inferieures a 1500 
Hz) ou I'influence de la tete sur le conditionnement de la matrice est assez 
moderee. 

Comme precedemment on peut constater qu'il n'existe pas d'ecartement 
angulaire permettant d'etre bien conditionne sur tout le spectre. Nous allons done 
etre inevitablement confronts a inverser la matrice de transfert electroacoustique 
dans des zones ou le conditionnement ne sera pas optimal. 
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3.2.4 La regularisation de la solution 


Kirkeby (8) propose alors de reconsiderer ce probleme d'inversion en termes 
de minimisation d'un critere quadratique. Dans ce but, on fait apparaitre un critere 
d'erreur que Ton va definir de la maniere suivante : 



Figure 17 : Modelisation du probleme en termes de minimisation de I’erreur quadratique 


Sur la figure 17, f est la frequence, T la periode d'echantillonnage, H le filtre 
optimal a determiner et C la matrice de transfert electroacoustique. 


L'erreur s'exprime ainsi comme la difference entre les 2 quantites suivantes : 

• Le signal original (u) multiplie par la matrice de decodage transaural (H) puis 
par la matrice de transfert electroacoustique (C) 

• Le meme signal original auquel on applique un retard de modelisation 
correspondant au temps de propagation entre les enceintes et I'auditeur ainsi 
qu'au retard induit par le filtrage des 2 matrices. 


Le critere de minimisation est alors defini dans I'espace de fourrier de la 
maniere suivante : 


J(k) = e*(k). e(k) + p. v*(fc). v(k) (15) 

Le premier terme correspond a notre erreur quadratique a minimiser (* est 
I'operateur de transconjugaison) mais on a aussi un terme de regularisation 3 qui 
pondere le signal de commande des haut-parleurs. 


On peut ainsi determiner le filtre optimal au sens des moindres carres : 

H(k) = [C H (k).C(k) + /?./d] -1 .C H (k).e n 


( 16 ) 
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On aller plus loin et choisir un parametre de regularisation qui depend de la 
frequence. Ceci permet alors de ne regulariser que la partie mal conditionnee de 
notre systeme (par exemple : >7kHz et <100Hz en configuration stereo-dipole). 


4 La problematique des basses frequences 


GENESIS utilise couramment des configurations a 4 haut-parleurs et un 
caisson de basse lors de Installation de systemes transauraux chez ses clients. La 
premiere paire de haut-parleurs est en general ecartee de 10° et la deuxieme de 
70°. Les frequences en dessous de 80 Hz ne sont pas transauralisees mais 
directement envoyees au caisson de basses, puisque ces frequences ne sont pas 
localisables. 


De tels systemes ne sont pas optimaux en termes de conditionnement, mais 
I'implementation d'un transaural en milieu industriel impose en general certaines 
contraintes : 

- La geometrie du systeme et son encombrement doivent s'adapter a la salle 
d'ecoute. 

- II faut pouvoir placer un ecran pour videoprojecteur entre les hauts 
parleurs frontaux pour que I'auditeur puisse interagir avec un logiciel, ce 
qui impose d'avoir un ecartement angulaire suffisant (souvent 10° ou 
plus). 

- Dans certains cas la localisation de sources arriere doit etre specifiquement 
favorisee pour certains types de sons automobiles, ce qui fonctionne 
empiriquement en ecartant Tangle pour la paire dediee aux basses 
frequences. 

- La position d'ecoute est parfois figee a des endroits non ideaux en termes 
d'influence des modes en basses frequences (presence de noeuds marques 
ou d'asymetrie consequente entre les 2 oreilles). 


Le cumul de ces contraintes vient souvent limiter les performances theoriques 
des systemes, et notamment dans les basses frequences. Or de nombreux clients 
de GENESIS viennent du secteur automobile ou aeronautique, ou les sons a 
restituer ont un contenu tres important dans cette zone. Legalisation finale en 
timbre du systeme ne suffit pas toujours a donner pleinement satisfaction au client 
lors des ecoutes, et la question se pose aujourd'hui des moyens qui peuvent etre 
mis en oeuvre par GENESIS pour ameliorer le rendu des basses frequences en 
restitution transaurale. 
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5 Presentation des differentes methodes envisagees 
pour I'optimisation des basses frequences 

Afin de repondre aux objectifs du stage, les 3 pistes suivantes ont ete 
investiguees: 

• Implementation du systeme transaural dit OSD (Optimal Source Distribution) 
sur 6 haut-parleurs. L'idee d'un tel systeme est de faire fonctionner des paires 
de haut-parleurs dans une gamme frequentielle adaptee pour optimiser le 
conditionnement de la matrice de transfert electroacoustique (chaque paire 
est done disposee avec un ecartement angulaire adapte). Contrairement a la 
methode 4.1 utilisee classiquement a GENESIS, les basses frequences sont ici 
incluses dans le processus de transauralisation. 

• Implementation de la methode d’inversion de Kirkeby pour le calcul du filtre 
d’egalisation final du systeme. Cette methode combine les techniques 
connues d'optimisation au sens des moindres carres et de regularisation de 
Tikhonov. Le parametre de regularisation devrait permettre davantage de 
souplesse sur I'elaboration du gabarit du filtre d'egalisation que la fonction 
FIR2 de Matlab utilisee a ce jour. 

• Implementation d'une methode d'adaptation de la matrice de transfert 
electroacoustique a I'ensemble {salle + enceintes}. Cette methode, dite de 
« pre-egalisation a 2 micros », consiste a prendre en compte la reponse des 
enceintes et I'effet des modes de salle aux positions des oreilles de I'auditeur. 

5.1 Le systeme OSD (Optimal Source Distribution) 

5.1.1 Principe du systeme ideal : 

Comme nous I'avons vu precedemment, I'utilisation du stereo-dipole sur la 
totalite du spectre audible entraine la necessite d'avoir a faire un compromis entre 
la dynamique, la robustesse et le controle que Ton a sur le systeme transaural. 
Nelson (9) a propose un systeme qui resout ces problemes fondamentaux en 
utilisant une distribution optimale des sources. 
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Figure 18 : Principe du systeme OSD 


On considers une paire conceptuelle de monopoles dont Tangle varie de 
maniere continue en fonction de la frequence : 



Figure 19 : Schema du systeme de deux monopoles pour 2 points d'acoute 


On a alors : 


/ nn \ 

0 = 26 = 2arcsin( — — — ) 

\2kAr) 

Ou O est /'angle entre les deux monopoles, k le nombre d'onde, Ar I'ecartement des 
oreilles, n un entier que I'on prend egal a 1 . 


( 17 ) 
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Dans le cas general il a ete montre que la matrice de transfert 
electroacoustique normalisee s'ecrit : 


C = 


ge 


-jkAl 


_ge 


-jkAl 


La matrice inverse normalisee peut alors s'ecrire : 


(18) 


H = C -1 = 


\ — g2 e -2 jk ArsinO ' 


ge 


ge 


-jkAl 


-jkAl 


(19) 


Ce qui met en evidence que le conditionnement depend de la frequence et 
qu'il y a des variations de dynamique dues aux variations de la norme de H avec la 
frequence. 


En se plagant dans le cas ideal de I'OSD ou I'ecartement de la paire de 
monopole suit la frequence selon I'equation (17), on obtient 


H = 


— f 1 

i + g 2 Y-jg 



( 20 ) 


On s'apergoit que ||H|| = 1/V2 et /(C) = 1 pour toute la gamme de frequence. 
Par consequent, il n'y a plus de perte de dynamique et I'erreur dans le calcul du 
filtre inverse est petite. Le systeme a done un excellent controle sur les signaux 
reproduits en plus d'etre plus robuste. 


5.1.2 Le systeme OSD discretise : 

En pratique, une paire de monopole dont Tangle varie continument comme 
une fonction de la frequence n'est pas realisable. Cependant il est possible de creer 
un systeme base sur ce principe en discretisant I'angle des haut-parleurs comme 
sur la figure 20. 
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a pair of discretised variable span OSD loudspeakers 



Le calcul du conditionnement de la matrice de transfert electroacoustique en 
fonction de la frequence realise precedemment permet de definir les angles et les 
bandes de frequences associees optimaux. En effet, pour un angle donne, la bande 
frequentielle ou I'amplification est relativement faible et ou la matrice de transfert 
electroacoustique est bien conditionnee, est relativement large autour de la 
frequence optimale. Par consequent cela nous autorise a limiter le nombre de paires 
de haut-parleurs sans pour autant perdre en controle des performances et en 
robustesse. 


Nelson choisit une configuration a 3 paires de haut-parleurs (10) (9): 

• La paire aigue est ecartee de 3.1° par rapport a I'axe de symetrie du systeme 
et est utilisee pour toutes les frequences superieures a 3500Hz 

• La paire medium est ecartee de 16° par rapport au meme axe et est utilisee 
pour toutes les frequences comprises entre 350 et 3500 Hz. 

• La paire grave est ecartee de 90° par rapport a I'axe et est utilisee pour 
toutes les frequences inferieures a 350Hz. 


Une representation du systeme est donnee en figure 21. 
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high 



Le decoupage a ete choisi de maniere a avoir un conditionnement optimal. En 
utilisant la figure du conditionnement d'un stereo-dipole (22) en fonction de la 
frequence, on obtient le schema de decoupage suivant. Chaque trait vertical 
represente le domaine frequentiel de travail de chaque paire de haut -parleurs. 



transducer span ( # ) 


Figure 22 : Domaines de frequences bien conditionnes pour les postions angulaires choisies 
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5.2 L'inversion de Kirkeby 


Lorsqu'un decodeur transaural est calcule, une phase d'egalisation finale est 
mise en place pour corriger le detimbrage du systeme. Un tel detimbrage est le 
resultat de I'effet de salle, du conditionnement non parfait a certaines frequences, 
et de I'approximation des HRTF par un modele de tete spherique. 


L'inversion de Kirkeby est une methode tres rapide de calcul de filtre inverse 
pour les systemes de reproduction audio. La methode combine I'optimisation au 
sens des moindres carres et la regularisation (11). On calcule un filtre FIR (finite 
impulse response) causal en optimisant la performance des filtres pour un grand 
nombre de frequences discretes. Par consequent, cette methode est utile seulement 
quand on a la possibility d'avoir des filtres inverses relativement long, ce qui est 
notre cas. 


On peut modeliser le probleme de deconvolution par le schema suivant : 

Target matrix 
x 


Recorded 

signals 

u(z) 


z 


Mz ) 




“7 


m) 


. Source 
input signals 

v(z) 


Matrix of 


Desired 

signals 

d(z) 



* Matrix of plant 

filters frarKfer fnnrfim 


Reproduced 
signals 


inverse filters transfer function w(z) 

Figure 23 : Schema du probleme de deconvolution multi-canal discretise en temps 


On se place dans un cas ou Ton cherche a recreer un signal d(z) a partir d'un 
enregistrement u(z). Pour realiser cela nous envoyons un signal v(z) a notre 
systeme de sortie. La matrice C(z), la matrice de transfert electroacoustique, et la 
matrice A(z) sont des donnees et I'objectif de la methode est de determiner H(z) en 
minimisant I'erreur e(z). 
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Pour determiner H(z), on definit une fonction cout, respectivement somme 
d'un terme d'erreur de performance (qui mesure la qualite de reproduction du 
signal) et d'un terme d'effort de penalite en suivant le principe de la regularisation 
precedemment decrit au paragraphe 3.2.4 : 


7(z) = e*(z). e(z) + fl.v*(z).v(z) 

* est I'operateur de transconjugaison 


( 21 ) 


Le parametre/? est un parametre de regularisation qui determine quel poids 
on attribue au terme d'effort. 


En faisant dependre ce parametre de la frequence on peut choisir 
d'appliquer une regularisation differente par bande de frequence et done limiter les 
problemes dus au mauvais conditionnement de la matrice de transfert electro- 
acoustique dans certaines bandes de frequences. 

5.3 La methode dite de pre-egalisation a 2 micros 

Comme mentionne au prealable, un modele spherico-ellipsoidal de tete est 
utilise pour determiner la matrice de transfert electroacoustique. Ce modele etant le 
resultat d'un calcul analytique effectue en champ libre, les HRTF qui en resultent 
sont totalement non bruitees, ce qui permet de s'affranchir du probleme d'echo 
present dans les reponses impulsionnelles lorsque ces HRTF sont directement 
mesurees (les salles d'ecoute cibles n'etant jamais anechoiques). 


Une egalisation finale en timbre du systeme permet done au systeme de 
s'adapter a la salle dans laquelle il se trouve. C'est notamment par cette egalisation 
que Ton va corriger (ou du moins attenuer) I'effet des modes de salle en basses 
frequences. Ces modes etant parfois tres prononces, de grandes excursions sont 
demandees au filtre final d'egalisation qui ne pourra les corriger que partiellement 
etant donne sa taille maximum autorisee (8192 points, pour des raisons de 
puissance CPU et d'effet « trainant » du filtre a I'ecoute). 


Afin de limiter ces excursions, I'idee de la methode proposee est de corriger la 
matrice de transfert electroacoustique theorique afin d'y inclure le detimbrage du a 
I'effet de salle sans y introduire d'echo. Ce faisant, le calcul de la matrice inverse 
aura deja pris en compte I'effet de salle, et le systeme devrait etre moins detimbre 
au moment de Legalisation finale. Celle-ci devrait alors etre plus souple, et I'effet 
des modes plus facilement corrige. Le principe de la methode est done le suivant : 


1/ On mesure avec un microphone AKG les reponses impulsionnelles de 
chacun des haut-parleurs aux deux points de I'espace correspondant aux oreilles de 
I'auditeur. 
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2/ On effectue un lissage en l/6eme d'octave du gabarit frequentiel des 
reponses (sinon les corrections hautes frequences seraient vite fausses des lors que 
I'auditeur s'ecarterait un peu de la position ideale). 

3/ On calcul des filtres a phase lineaire a partir de chacun des gabarits 
frequentiels lissees. 

4/On vient convoluer les reponses impulsionnelles des filtres ainsi mesures 
avec chacun des transferts electroacoustiques modelises dans la matrice initiale. 

5/ On effectue le calcul du decodeur transaural comme avant, mais a partir de 
cette nouvelle matrice de transferts electroacoustiques. 

6 Resultats obtenus 


6.1.1 Apport de I'OSD 

II s'agit dans cette partie de comparer les resultats obtenus avec configuration 
OSD (cf annexe) a ceux qui correspondent a une configuration classique en 4.1 
employee chez GENESIS, a savoir : 

• Une paire avec ecartement angulaire de 10° pour restituer le son sur la plage 
3kHz-20kHz 

• Une paire avec ecartement angulaire de 70° pour restituer le son sur la plage 
70Hz-3kHz 

• Un caisson de basses pour restituer les frequences en dessous de 70 Hz 


6. 1.1.1 Apport de I'OSD sur le timbre final 

Le rendu du timbre est une exigence importante dans les systemes 
transauraux puisque c'est ce qui rend I'ecoute fidele a I'enregistrement. L'utilisation 
de I'OSD permet d'obtenir de tres bons resultats sur la reponse en frequence finale 
du systeme. Comme on peut le voir sur les figures ci-dessous (25 et 26), la 
configuration a 4 haut-parleurs presente des variations de +/-4dB dans les HF. 
Ponctuellement il peut y avoir jusqu'a 3 dB d'ecart entre les deux oreilles et on note 
des pertes importantes en basses frequences. L'utilisation de I'OSD permet : 

• Une reduction des variations a +/-2dB. 

• Un gommage des ecarts de niveau entre les 2 oreilles. 

• Un rendu des basses frequences fidele jusqu'a 30dB. 
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Reponse en frequence finale 4HP 



Figure 24 : Reponse en frequence du systeme transaural a 4 haut-parleurs et un caisson 


Reponse en frequence finale 6HP 



Figure 25 : Reponse en frequence du systeme transaural a 6 haut-parleurs 


6. 1.1.2 Rejection 

La rejection du systeme caracterise I'aptitude a annuler les chemins croises 
sur toute la gamme de frequence. Elle est done d'une importance primordiale pour 
la spatialisation du son puisqu'il faudrait idealement que chaque oreille ne regoive 
qu'un seul canal de I'enregistrement binaural. 
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Pour realiser une mesure de rejection nous jouons un bruit rose a gauche puis 
a droite et nous comparons I'ecart de niveau entre les deux oreilles a chaque fois 
pour obtenir respectivement la rejection a droite et la rejection a gauche. 


En basses frequences la configuration a 4HP ne transauralise pas. On peut 
toutefois observer la presence d'un mode a 40Hz qui implique des differences 
importantes entre les deux oreilles (6dB d'ecart). L'OSD permet d'apporter la 
transauralisation avec une rejection moyenne de 5dB. On retrouve le mode a 40Hz 
qui a ete amplifie par la transauralisation. Nous verrons par la suite qu'une pre- 
egalisation du systeme permet de corriger un peu cet effet d'amplification. 


A 3000Hz, nous observons que la configuration a 4HP donne une rejection 
positive. Ceci est probablement le resultat d'un probleme de conditionnement a 
cette frequence, ou d'un effet d'une reflexion parasite. En effet une position 
angulaire de 10° des haut-parleurs frontaux ne permet d'avoir un bon 
conditionnement dans cette zone de frequence, et la salle est legerement traitee 
acoustiquement mais loin d'etre anechoique. L'OSD comme nous pouvons le 
remarquer n'a pas ce probleme puisque la disposition des haut-parleurs cherche a 
optimiser le conditionnement de la matrice de transfert electroacoustique pour 
faciliter son inversion. 


En hautes et moyennes frequences on obtient des rejections de I'ordre de 10 
dB avec I'OSD alors que la configuration a 4 haut-parleurs presente une rejection 
moyenne inferieure a 5dB. 


Rejection Oreille Gauche 4HP Rejection Oreille Droite 4HP 



Figure 26 : Rejection de I'oreille gauche et de I’oreille droite dans le systeme transaural a 4 haut-parleurs et un 
caisson en fonction de la frequence 
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Rejection Oreille Gauche 6HP Rejection Oreille Droite 6HP 




Figure 27 : Rejection de l'oreille gauche et de I’oreille droite dans le systeme transaural a 6 haut-parleurs en 
fonction de la frequence 


L'ecoute comparative des deux systemes par les employes de GENESIS donne 
I'OSD meilleur en spatialisation et en timbre comme attendu. II faudrait neanmoins 
realiser une etude psycho-acoustique complete pour tirer de plus amples 
conclusions. 


6.1.2 Apport de Kirkeby dans I'inversion du filtre d'egalisation 

Comme nous I'avons vu theoriquement, la methode d'inversion de Kirkeby 
permet d'avoir plus de controle sur I'inversion du filtre final d'egalisation et done 
d'ameliorer le timbre final du systeme transaural. Cette methode permet via la 
regularisation par bande de frequence de stabiliser I'inversion des basses 
frequences et d'eviter I'apport d'energie en hautes frequences. 


Sur la figure 28 nous observons en vert le filtre d'egalisation obtenu lorsqu'on 
n'utilise aucune regularisation sur tout le spectre. La reponse du filtre est assez 
instable en basses frequences, on observe des oscillations plus ou moins marquees 
jusqu'a plus de 3000Hz. De plus la moyenne de ses oscillations ne suit pas 
exactement le filtre cible attendu et observe des ecarts de plus de 2dB (et meme 
ponctuellement des ecarts de 8dB) ce qui est inacceptable compte tenu de la 
precision recherchee. 
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Figure 28 : Calcul du filtre d'egalisation final sans appliquer de regularisation a la methode d'inversion de Kirkeby 


Calcul du filtre d egalisation par la methode de Kirkeby avec une bonne regularisation 



Frequences (Hz) 


Figure 29 : Calcul du filtre d'egalisation final en choisissant bien les parametres de regularisation de I'inversion de 

Kirkeby 
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En choisissant bien les parametres de regularisation sur tout le spectre, nous 
parvenons a stabiliser le filtre et supprimer les oscillations. De plus le filtre presente 
des ecarts moins importants par rapport au filtre cible. 


6.1.3 Apport de la pre-egalisation 

Comme nous I'avons vu, la pre-egalisation permet d'introduire les effets lies a 
la salle ainsi qu'au systeme de restitution dans la matrice de transfert electro- 
acoustique avant de I'inverser en utilisant les gabarits frequentiels des reponses 
mesurees aux deux positions des oreilles par des micros AKG avec le modele 
spherico-ellipsoidale. 


Nous regardons la reponse du systeme transaural avant egalisation pour 
observer I'apport de la dite pre-egalisation : 


Reponse du systeme transaural avant egalisation 



Figure 30 : Reponse en frequence du systeme transaural a 6 haut-parleurs avant egalisation 


On observe des ecarts importants allant jusqu'a 8dB entre les deux positions 
d'ecoute. On remarque aussi egalement la presence du mode a 40Hz 
precedemment decrit. De plus on note la difficulty qu'a le systeme a restituer les 
tres basses frequences (le caisson ne descend naturellement qu'a 40 Hz). 
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Reponse du systeme transaural corrige avant egalisation 



Figure 31 : Reponse en frequence du systeme transaural a 6 haut-parleurs pre-egalise avant egalisation 

En realisant la pre-egalisation nous obtenons la figure 31. On note que les 
differences de niveau interaurales sont corrigees sur la quasi-totalite du spectre. La 
reponse globale du systeme est plus plate. L'ensemble de ces parametres nous met 
dans de bonnes conditions pour realiser un filtre d'egalisation efficace par la 
methode de kirkeby. 


Line fois cette egalisation finale realisees nous obtenons la reponse en 
frequence ci-dessous : 


Reponse en frequence finale 6HP 



Figure 32 : Reponse en frequence du systeme transaural a 6 haut-parleurs 


Yannick Busawon 
Avril 2013 
Page 39 / 45 


jqenesis 

H “all dimensions of Sound 


GENESIS 

Restitution des Basses Frequences et optimisation 
dans les systemes Transauraux 3D 


Ici, la reponse du systeme est tres homogene et les differences entre oreilles 
sont gommees. Alors que si nous observons la figure 34 qui represente la reponse 
en frequence finale sans pre-egalisation, nous notons que il reste des bandes de 
frequences mal egalisee et les basses frequences s'effondrent a partir de 50Hz alors 
que nous gardions une reponse plate jusqu'a 30Hz dans le cas avec pre-egalisation. 


Reponse en frequence finale 6HP sans pre-egalisation 



Figure 33 : Reponse en frequence du systeme transaural a 6 haut-parleurs qui n’a pas ete pre-egalise 


En termes de rejection, la pre-egalisation permet d'attenuer le mode a 40Hz 
sans pour autant parvenir a le corriger totalement. On observe aussi une rejection 
meilleure sur I'ensemble du spectre. 


Rejection Oreille Gauche 6HP sans pre-egalisation Rejection Oreille Droite 6HP sans pre-egalisation 




Figure 34 : Rejection de I'oreille gauche et de I’oreille droite dans le systeme transaural a 6 haut-parleurs non pre- 
egalise en fonction de la frequence 
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Rejection Oreille Gauche 6HP 



Figure 35 : Rejection de I'oreille gauche et de I’oreille 
fonction de 


Rejection Oreille Droite 6HP 



Frequences (Hz) 

droite dans le systeme transaural a 6 haut-parleurs en 
la frequence 
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7 Conclusion et perspectives 


L'objet du stage etait de realiser un systeme permettant de transauraliser 
toute la gamme audible tout en ameliorant la qualite du rendu en timbre et en 
rejection. Pour cela nous avons etudie et realise un systeme transaural OSD 
(optimal source distribution) a 6 haut-parleurs, avec un filtre d'egalisation calcule 
par la methode de Kirkeby et ce que nous avons appele la pre-egalisation du 
systeme. 

L'OSD nous a permis de realiser un transaural sur tout le spectre car le choix 
des angles entre haut-parleurs nous permet de minimiser les defauts d'inversion 
dus au mauvais conditionnement de la matrice de transfert acoustique. En 
s'affranchissant de ces problemes nous obtenons alors une rejection sur I'ensemble 
du spectre de I'ordre de lOdB la ou nous obtenions au mieux 5dB avec I'ancien 
systeme. Le rendu perceptif est tout a fait satisfaisant en termes de spatialisation 
notamment grace a cela. 

L'inversion de Kirkeby nous a permis de realiser une meilleure egalisation du 
systeme. Le timbre est particulierement important dans le cahier des charges des 
clients de GENESIS. Nous avions des ecarts de +/-5dB avec le systeme a 4 haut- 
parleurs et nous sommes passe a +/-2dB ce qui est entierement satisfaisant. En 
plus de cela nous corrigeons les basses frequences plus efficacement jusqu'a 30dB. 


La pre-egalisation quant a elle nous a permis de corriger les defauts de la 
salle et de notre systeme qui engendraient des ecarts notables entre oreilles. A cela 
il faut ajouter que cela nous permet de faciliter le calcul du filtre d'egalisation car les 
reponses en frequences sont pre-egalisees. 


Le systeme transaural est aujourd'hui bien meilleur a la fois timbre et en 
rejection mais on peut imaginer d'autres ameliorations venant parfaire le resultat : 

- On peut utiliser des caissons de basse supplementaires pour corriger les 
modes de salle et obtenir des resultats meilleurs en tres basses 
frequences. 

- On peut imaginer ^implementation d'un filtre d'egalisation en temps en reel 
avec un HR (infinite impulse response) 
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Annexe 



Figure 36 : Configuration OSD realisee 


La photo represente la configuration mise en place pendant mon stage. On peut noter 
qu'il existe un traitement acoustique des parois avec de la mousse mais cela n'est pas 
efficace pour les basses frequences et on observe de forts modes de salle. 


Les enceintes sont ecartees de 3°, 16° et 90° par rapport a I'axe central. 


